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Summary. We considered the Edwards- Wilkinson model on a small-world net- 
work. We studied the finite-size behavior of the surface width by performing exact 
numerical diagonalization for the underlying coupling matrix. We found that the 
spectrum exhibits a gap or a pseudo-gap, which is responsible for a finite width 
in the thermodynamic limit for an arbitrarily weak but nonzero magnitude of the 
random interactions. 



1.1 Introduction 

Since the introduction of small- world networks [1,2], a number of well-known 
models have been studied where the original short-range interaction topol- 
ogy is extended to include a finite number of possibly long-range links per 
"site" [3]. The common observation is that these systems can undergo a 
phase transition, even when the random links are added (or rewired) to a 
one-dimensional original substrate [3-8] . The nature of the transition resem- 
bles that of a mean- field one [6-8] . 

Among the first applications of small-world networks was to study syn- 
chronization in dynamical systems [2,9] such as the Kuramoto oscillators [10]. 
The need for autonomous synchronization for a system with a large number of 
"agents" (processing elements) also naturally emerges in large-scale parallel 
discrete-event simulation (PDES) schemes [11] for systems with short-range 
interactions and asynchronous dynamics [12,13]. Frequent and necessary "lo- 
cal" communications between processing elements (PEs) to ensure the asyn- 
chronous causal dynamics of the underlying system will eventually lead to a 
diverging spread of the progress of the individual PEs [13, 14]. This property 
can seriously hinder efficient data collection for such simulation schemes. An 
alternative to possibly costly and frequent global synchronizations is to ex- 
tend the required short-range communication topology to include "weak" ran- 
dom links [14]. Weak in this context refers to the relative timescale of actually 
using the random connections for synchronization. By directly "simulating 
the simulations" and using simple coarse-graining arguments, it was demon- 
strated [13] that the progress of the PEs with only local synchronization 
exhibits "kinetic roughening" governed by the Kardar-Parisi-Zhang (KPZ) 
equation. With random links added (a finite number per PE) and invoked at 



2 B. Kozma, G. Korniss 



an arbitrarily small (but non- vanishing) rate, however, the PEs progress in a 
near- uniform fashion [14]. 

Here we focus on how critical fluctuations (originally present in the steady 
state of a one-dimensional system) are suppressed when the interaction topol- 
ogy is extended to include weak interactions facilitated by random links. To 
this end we study the Edwards- Wilkinson (EW) linear stochastic growth 
equation on a "substrate" with small-world-like topology. This model is also 
closely related to phase ordering and synchronization among coupled oscil- 
lators in the presence of noise [8] and to the XY-model on a small-world 
network [6]. We consider the equation 



where hi is the surface height, T]i(t) is a delta-correlated Gaussian noise with 
variance 2 (without loss of generality), and we have dropped the i-dependence 
from the argument of hi for brevity. The matrix represents the (quenched) 
random links on top of a one-dimensional lattice of length N (even for sim- 
plicity) with periodic boundary conditions, i.e., Jy=l if a random link is 
present and zero otherwise. The parameter p is the strength of the interac- 
tion through the random links. Our construction of the random links is such 
that each site has exactly one random link. More specifically, pairs of sites are 
selected at random and once they are chosen, they cannot be selected again. 
This somewhat constrained construction of the random network originates 
from an application to scalable PDES synchronization schemes [14], where 
fluctuations in the individual connectivity of the PEs are to be avoided. 

For a given realization of this small-world network the average surface 
width characterizing the roughness is defined as 



where h=(l/N)'^2^ r _ 1 hi is the mean height and (. . .) denotes an ensemble 
average over the noise in Eq. (1.1). For p~0 in Eq. (1.1), we recover the one- 
dimensional EW model where the steady-state width diverges as (w%)=N/12. 

One may wonder how the system would behave if the same total number 
of links as in the above construction of a small- world network (i.e., N/2) were 
used to connect each site with the one located at the "maximum" possible 
distance of away from it (N/2 on a ring with periodic conditions). Elementary 
calculations show that (w%) ~ iV/24 for large N, i.e., the width would diverge 
as for a one-dimensional system of size N/2. Indeed, one can realize that 
such regularly patterned long-range links make the original system equivalent 
to a 2x(N/2) system with only nearest-neighbor interactions and shifted 
periodic boundary conditions. More generally, one can show that, if every 



N 



d t hi = -(2hi - h i+ i - hi-x) - p^2j i:j (hi - hj) + r]i(t) , (1.1) 



i=i 




(1.2) 
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site is connected to a finite number of others and the length of those links 
can only assume a finite set of long-range values (all scaling with N), the 
width will still diverge in the same fashion as for the one-dimensional case. 



1.2 Roughness and the Density of States 

We study the finite-size effects of the width of the surface and also the un- 
derlying spectrum (density of states) of the associated random matrix which 
governs the steady-state height fluctuations. Exploiting that the noise in 
Eq. (1.2) is Gaussian, the steady-state width for a single realization of the 
random network can be expressed as 

where Xk are the eigenvalues of the real symmetric coupling matrix 

Jij = {(2 + p)Sij - S i+ i j - Si-i j} - pJij (1.4) 

as can be read off from Eq. (1.1). Note that we have exploited our specific 
choice of Jy, resulting in Jij = ^ f° r au Also, note that since Eq. (1.2) 

contains the height fluctuations measured from the mean, the eigenvalue 
Ao=0, corresponding to the uniform eigenvector (zero-mode) of r^, does not 
appear in Eq. (1.3). In the limit of N^>oo and assuming that the distribution 
of the eigenvalues of becomes self-averaging, the disorder-averaged width 
can be be written as 



N-l 



1^1_ 



N^oo f p {X)dX ^ 
X 



where [. . .] stands for averaging over the random-link disorder and p(X) de- 
notes the density of eigenvalues of F^. The behavior of p(X) as A goes to 
zero determines whether the width remains finite or diverges in the thermo- 
dynamic limit. In the pure one-dimensional case, p(A) actually diverges as 
l/(2ny/X). If, however, p(X) exhibits a gap or approaches zero fast enough, 
[(u> 2 )tv] will be finite. In the context of diffusion on a small-world network, 
it was found that the density of states exhibits a pseudo-gap (vanishes ex- 
ponentially fast) [15]. The construction of the small- world graph in Ref. [15] 
allowed for the existence of arbitrarily long "pure" chain-segments of the net- 
work with exponentially small probabilities. These small, but non-vanishing, 
probabilities were responsible for the pseudo-gap [16]. In our specific con- 
struction of the network, where each site has exactly one random link, the 
above argument does not apply and a true gap may develop. Further, the 
coupling matrix [Eq. (1.4)] has a realization-independent "mass" term. This 
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Fig. 1.1. (a) Disorder-averaged surface width [(lOjv)] as a function of the system size 
N for p = 1, 1CT 1 , 10~ 2 , 10~ 3 , 10~ 4 , and 10~ 5 (from bottom to top, respectively, 
with symbols), obtained by using 1000 independent realizations of the random 
network (except for iV=2000 where only 100 realizations were generated). The solid 
lines are the analytic forms for the width in the simple "massive" approximation. 
The straight line corresponds to the "rough" p—0 case, (b) Cumulative (integrated) 
density of states for p — 10 _1 , 10~ 2 , 10~ 3 , and 10~ 4 (from right to left, respectively, 
with bold solid lines) based on 1000 realizations of the random network for iV=1000. 
The thin lines correspond to the analytic form for the infinite system-size "massive" 
approximation. These curves are not distinguishable on normal scales, except for 
the p=W~ 1 case. The inset shows the same on log-log scales to magnify the region 
near small eigenvalues. We also plot the analytic form for the infinite-system p=0 
case (asymptotically a straight line [~ (l/7r)vA)] for small A. 

property would actually allow for a perturbation expansion for small but 
non-zero values of p with the term —pJij being the perturbation. 

We performed exact numerical diagonalization of the coupling matrix 
Eq. (1.4) using standard numerical routines [17], and calculated the steady- 
state width as a function of the system size for various values of p. The results 
are summarized in Fig. 1 (a) . We also plotted the analytic form of the width for 
the simple "massive" coupling matrix, the expression in brackets in Eq. (1.4), 
as the zeroth-order approximation in a pcrturbative approach. It appears 
that for small values of p, the numerically computed (and disorder-averaged) 
width and this simplest approximation yield the same asymptotic finite-size 
effects. Fig. 1(b) shows the cumulative eigenvalue distribution J p(X')d\' 
for TV^IOOO for various p values. Whether the spectrum exhibits a true or 
a pseudo-gap (due to exponentially small likely eigenvalues), cannot be de- 
termined by numerics. It is appears, however, that the numerically observed 
"gap" asymptotically scales linearly with small values of p. 
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I. 3 Conclusions 

We carried out exact numerical diagonalization for the coupling matrix rep- 
resenting EW growth on a small-world network. In our construction each site 
had one random link, i.e., no fluctuations were allowed in the connectivity. We 
found that the surface width saturates for all nonzero values of the amplitude 
of the random coupling as a result of the gap or pseudo-gap in the under- 
lying spectrum. We should also note the similarity between the relaxation 
properties of our model and that of the one-dimensional Ising-like systems 
with ("annealed") random spin-exchange process [18] . This long-range pro- 
cess creates a mean-field-like environment, in which ordering is possible with 
other suitable chosen local processes present [18]. 
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